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^ Abstract 

<N 

Either by solving the ghost propagator DSE or through a one-loop computation in the 
^ (— I RGZ (Refined Gribov-Zwanziger) formalism, we show that a non-trivial ghost-gluon vertex 

pi I is anyhow required to obtain a ghost propagator prediction compatible with the available 

' corresponding lattice data in the SU(3) case. For the necessary gluon propagator input, we 

present RGZ tree level fits which account well for the gluon lattice data. Interestingly, this 
^ propagator can be rewritten in terms of a running gluon mass. A comparison of both DSE 

and RGZ results for the ghost propagator is furthermore provided. We also briefly discuss 
the connection between the RGZ and the OPE d = 2 gluon condensate. 
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O 1 Introduction 

(N 

. • As the Green functions of a general quantum field theory contain in principle all information, 

considerable research efforts is being devoted to construct as reliable as possible estimates of 
^ these Green functions. Amongst the various possible ways to compute the Green functions, 

^ a powerful method is the study of the Dyson-Schwinger equations (DSEs), which are nothing 

else than the quantum version of the equations of motions and are thus derivable from the 
partition function itself [1-5]. The appealing it may look to actually solve the DSEs, it must be 
mentioned that they inevitably lead to an infinite number of coupled equations and consequently 
assumptions/truncations must be made to obtain any kind of results. For SU(iV) Yang-Mills 
gauge systems, usually quantized in the Landau gauge because of the special (renormalization) 
properties of this gauge, the DSEs for the propagators themselves are solved while the input 
vertices are either taken to be tree level like or modeled using information from e.g. lattice 
simulations of the non-perturbative vertices when available [4, 6-14] . 

The Landau gauge is also perfectly suited for a numerical lattice approach, allowing for 
explicit comparisons between continuum results and their lattice counterparts (for a recent 
review see [15]), which are a priori non-perturbative and exact, apart from the usual difficulties 
of correctly handling discretization artifacts [16-20]. This means certain assumptions can be 
tested, for example if it is allowed to invoke for certain vertices the tree level value. The here 
presented work fits herein. Specifically, we will discuss the solution of the ghost propagator DSE 
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(GPDSE), using a fine modeling of the input gluon propagator from precise large volume lattice 
data. The variable ingredient of the DSE will be the ghost-gluon vertex. 

In a previous paper [21], some of us already discussed how to model the ghost-gluon vertex 
inspired by an Operator Product Expansion (OPE) analysis with a d = 2 condensate, allowing 
for a reasonable description of the available vertex' lattice data. With the exception of some 
(semi-classical) analysis' of the vertex made in [13, 22] and the very recent d = 2 work of [23], 
other DSE works usually deploy a tree level ghost-gluon approximation. As advocated in several 
independent series of DSE papers, qualitative agreement between the gluon/ghost lattice data 
is found with what is now known as the massive gluon (decoupling) DSE solution (i.e., see [4, 8- 
10]). This qualitative agreement should be understood as that the coupling at the renormalization 
point is rather treated as a fitting parameter in order to account for the lattice data. Otherwise, 
mismatches between DSE and lattice results appear, as it is particularly manifest for the ghost 
dressing function [24]. It should be remembered that, after the DSEs have been appropriately 
renormalized, the renormalized coupling could also be taken from the lattice as an input, in 
which case clearly is not an additional fitting parameter. Although such a qualitative agreement 
may be the first concern, once this is settled one should also aim at agreement at the quantitative 
level, otherwise we cannot claim that the DSEs capture all relevant dynamics. 

In this paper, we set an important step in this program. We shall show that a tree level (bare) 
ghost-gluon vertex never leads to an acceptable ghost DSE estimate of the lattice ghost propa- 
gator. Such could already have been guessed from the results of studies in this approximation, 
in which case there was always a mismatch with the lattice ghost (sec e.g. [8,24,25]). Incor- 
porating a nontrivial vertex model, compatible with the lattice data, based on the preparatory 
work [21] does lead to a fine match. Simplified versions of this vertex model also keep leading to 
discrepancies in the lower momentum regime, as well do the analytical one-loop approximations 
to the ghost propagator using the Refined Gribov-Zwanzigcr formalism, again using the latter 
simplified models. This clearly illustrates the important role played by the ghost-gluon vertex 
to adequately describe the Landau gauge ghost propagator within the DSE approach. 

In what follows, we shall first describe our setup and used input, before coming to results 
and our conclusion. 

2 A precise modeling of the input gluon propagator 

The main motivation of this paper is the analysis of the ghost dressing results obtained from the 
ghost DSE, particularly in connection with the input for the ghost-gluon vertex. The first step 
should be to provide ourselves with a proper and reliable model for the gluon propagator, which 
accounts for the corresponding "ab initio" QCD lattice data. As will be discussed below, this 
will permit us to solve the ghost DSE (in isolation from the gluon DSE) with the ghost-gluon 
vertex, a main target of this study, as the only unknown ingredient, and finally compare with 
available ghost lattice data. 

2.1 The RGZ gluon propagator and a running gluon mass 

Since the key work of Gribov [26], it has become well-appreciated that the Landau gauge does 
not allow to pick out a single gauge field per orbit: multiple gauge equivalent fields all fulfilling 
d^Af^ = can exist. As a consequence, the usual Faddeev-Popov procedure becomes in principle 
inadequate at the non-perturbative level. The Gribov copy problem is indeed not expected to 
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play a role in the ultraviolet where perturbation theories applies due to asymptotic freedom. 
However, the infrared dynamics is severely changed. 

In the continuum, the up-to-date best worked out (partial) cure to the Gribov problem 
— at least to our knowledge — is based on the original Gribov-Zwanziger approach: when the 
functional 7^[yl] = J di^xA^ is minimized per gauge orbit, a necessary condition for an cxtrcmum 
is the transversality condition dfj,A^ = 0, i.e. the Landau gauge. A local minimum is approached 
when the Faddeev- Popov operator, —df^D^, is positive The latter requirement is closely related 
to the occurrence of Gribov copies linked with zero modes of the Faddeev-Popov operator, as it 
is immediate to recognize that the existence of a mode oj with —d^D^u = leads to the existence 
of an (infinitemisal) gauge copy due to [A^ -|- eD^w] = d^j^A^. Therefore, as a first (and to date 
the only practically worked out) way to at least kill off this subset of infinitesimally related gauge 
copies in the continuum world, Gribov and Zwanziger implemented the restriction of the gauge 
path integral to the so-called Gribov region, where the Faddeev-Popov operator is positive^. 
Let us refer to the original works [26, 29] and the recent overview [30] for the underlying details. 
For the purpose of this paper, it is sufficient to know that this restriction can be brought into 
the working form of a new partition function, entailing the presence of a non-perturbative mass 
(Gribov) parameter 7^ oc A^^^. Formally, for 7^ = 0, the conventional Faddeev-Popov action 
is recovered, so it becomes evident that the Gribov copies' effect, as carried by 7^, is a pure 
non-perturbative effect, and the presence of the "soft" mass parameter 7^ is only infiuencing 
the infrared structure of the theory. For example, the predicted tree level gluon propagator was 

Due to the apparent mismatch of this propagator with state-of-the-art lattice data of e.g. [31- 
37] , it became clear this was not the end of the story. In a series of papers [38-41] , it has been 
discussed that the Gribov-Zwanziger action dynamically improves itself by the condensation of 
dimension 2 condensates, including the elsewhere much studied gluon condensate {A^) {e.g., see 
[42-48]). The tree level gluon propagator of this Refined Gribov-Zwanziger (RGZ) scheme is 
provided by 

I - p4 + + ^2)^2 + ^4 ' (2) 

where the parameters and are related to the local condensates of the gluon and auxiliary 
fields intended to preserve locality of the RGZ action [39-41]. In particular, 

- . (3) 

where 5'^(^^)rgz is the dimension- two condensate of the local operator A"^, also appearing in 
the OPE analysis of lattice Green functions [42,49-51]. The renormalization prescription for 
this local operator followed in OPE or RGZ approaches may differ and, as will be discussed 

'^This statement makes sense in the Landau gauge in which case the Faddeev-Popov operator is Hermitian. 
^In the lattice approach to Landau gauge dynamics, the Gribov problem is also there, see e.g. [27] and [28] 

for a recent toy model analysis, and it is usually treated in a quite similar fashion: one restricts the gauge fields 
admissible to the Monte Carlo sampling to a subset belonging to the first Gribov horizon, i.e. one numerically 
searches for the "optimal" minima of the Landau gauge functional. 
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below, this would imply that the non-perturbative condensates are differently defined. The A 
parameter in Eq. (2) can be also related to the Gribov parameter 7^, 



= m^M^ + 2g'^Nc-i^ . (4) 

Eq. (2) has been successfully apphed to describe SU(3) and SU(2) lattice gluon propagators [52, 
53] . We shall therefore also use it here to model in an adequate fashion the MOM-renormalized 
gluon propagator. It is interesting to notice that, as also done in [54], the propagator (2) can 
be reinterpreted as 

in terms of what we can call an effective "running RGZ gluon mass"^, 
with 

and where \j} is the renormalization momentum at which the MOM prescription determines 
that 

.(^^) = 1 + =%S. (8) 
A* 

It should be noticed that the mass function in Eq. (6) is gifted (as can be also seen in Fig. 2) 
with the main feature of reaching a positive saturation point (frozen mass), at asymptotically 
low momenta, that shifts toward a negative correction when the momentum increases. 



ml for fc^ < 



This feature has also been — at least qualitatively — observed in the results from lattice simula- 
tions. In particular, the OPE power corrections to the perturbative prediction for the gluon 
propagator, that have been shown as unequivocally needed to describe the lattice data in 
refs. [42,49], appears to behave as Eq. (5) with Eq. (9) and dominated at intermediate mo- 
menta by the dimension-two gluon condensate of the same local operator. This is a qualitatively 
striking result, although, as will be seen below, the condensates in RGZ and OPE approaches 
might have been differently defined. 

Then, Eq. (5) is a good candidate to model the gluon propagator and will be confronted to 
SU(3) lattice data in the next subsection. Before this, a few words about how the condensates 



^For a discussion on the running gluon mass and its modelling see e.g. [13,55-57]. 
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are defined and can be related in both approaches might be here in order. On should remember 
that, in the OPE analysis, one would be left with 



where all the quantities are assumed to be renormalized at the subtraction point fj? and, at 
tree-level, one obtains [42] 



The local operator contains a quadratic UV divergence that should be properly subtracted 
(see, for instance, the discussion of ref. [58]), for instance, by introducing a normal product 
referred to the QCD perturbative vacuum. In particular, the non-perturbative condensate in 
Eq. (10) should be specifically understood such that, given a l.h.s. that is non-perturbatively 
computed {e.g., on a lattice) and the Wilson coefficients cq and C2 computed in perturbation 
at a certain order and for a particular renormalization scheme, as the product of C2 and the 
condensate should be free of divergencies^, {A'^) is unambiguously defined. 

In the underlying effective potential approach to (A^)^^,^, the divergences of are consis- 
tently been taken care off using an order by order multiplicative renormalization [41, 45, 46]. As 
discussed in [45], this approach allows to investigate non-perturbative UV contributions to {A"^). 
If (^4^) is not a purely infrared-physics-dominated condensate, part of it escapes the usual OPE 
dictionary. This observation reinforces the fact that (A^)^^^ and {A^)^^^^ are not necessarily 
the same objects^. 

Since the exact computation of the effective potential is out of the question, as is thus the 
exact determination of the condensate(s), one can try to use the lattice data as a guide towards 
the value of e.g. (A^)^^^. We can also expand the RGZ gluon propagator from Eq. (5) at large 
momenta in powers oi l/k^ and, after comparing with Eq. (10) at tree-level, we will be left with 



This last result, grounded in the tree-level RGZ and OPE expansions, has to be taken of course 
as a rough approximation, but it already gives a clear indication that both condensates are 
differently defined, and do play a different role. In the OPE approach {A^)qpe serves to incorpo- 
rate non-perturbative infrared effects in the high(er) mom,entum, regime, while (^^)rgz and the 
other RGZ condensates arc relevant for the low(er) momentum, regime. As we will mainly focus 
on the tree level RGZ propagator here, we are clearly omitting any log contributions which are 
known to be important once the momentum gets larger^. So, we should not expect a tree level 
RGZ propagator to fit the whole momentum region. We shall come back to this subtle issue 
shortly. 

*In the Shifman-Vainshtein-Zakharov sum-rules approach [59, 60] , the UV divergences for the local operator 
and the IR ones coming from the Wilson coefficient cancel each other. 

^As there was some recent controversial work about the notion of condensates (sec [61,62] and references 
therein), we can only mention the condensates in the current work have to be understood as mass-scale parameters 
respectively related to the gluon propagator OPE expansion or the RGZ effective potential. Trying to properly 
settle the connection of both condensates (and check this from lattice results) remains as an interesting open 
question. 

®At low momentum, the mass scales will shield the Landau pole, allowing for a "freezing" of the logs. 
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10.89 


14 



Table 1: Lattice setup. The number of configurations refers to those used to compute the gluon propa- 
gator. For the ghost propagator on the 80"^ lattice only 11 configurations were considered. The physical 
scale was set using ro = 0.5 fm. We refer to [36] for more details. 



2.2 The SU(3) lattice gluon propagator 

In the present work, we will borrow the lattice data simulated at very large lattices (lattice 
lengths ranging from around 11 to 16 fm) by the Berlin-Moscow- Adelaide group [36]. Of all 
available volumes, we will consider only those where data is available both for the gluon and 
ghost propagators. 

The lattice QCD simulations were performed on a finite 4D torus using the Wilson action. 

For pure Yang-Mills simulations one can take the lightest glueball mass ~ 1.7 GcV as a typical 
hadronic scale, the corresponding length scale being L ~ Q.l fm. In order to investigate the 
infrared properties of the theory one should use sufficiently large volumes to avoid the problem 
of finite size effects and a resolution which accommodates the typical length scales for pure 
Yang-Mills theory. In practice, we have to compromise between using a large physical volume 
and a minimal length defined for the lattice. For the data sets considered here, the physical 
volume is sufficiently large to expect essentially no finite volume effects. However, the lattice 
spacing a 0.17 fm is about twice the typical hadronic length scale and one can expect finite size 
corrections due to the relative large a. For the gluon propagator the analysis of the combined 
effects of finite volume and finite lattice spacing is described in [20]. The ensembles considered 
in this work are described in Tab. 1. 

If we apply a MOM renormalization prescription at the subtraction point /j? , we get 



where £)Latt|-^2-| jg ^^iq bare lattice propagator and the renormalization constant is non- 

perturbatively defined as 

r)MOMc 2x1 _ _ jL (14) 

This last renormalization prescription is required (13) to be fulfilled at the subtraction point 
H = 3.6 GeV. 

As the main piu-pose of this letter is to investigate the ghost propagator Dyson-Schwinger 
equation, we need an as good as possible fit to the lattice gluon data over a so wide as possible 
momentum regime. The free RGZ-model inspired formula in Eq. (5) can do this job. This is 
achieved by introducing a global rescaling factor Z as an additional fitting parameter for the 
tree-level RGZ propagator given by Eq. (2), 

r>( ■2\ - 7 p^ + M^ 
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Figure 1: Renormalized gluon propagator D{p^) (left) and gluon dressing function p^iD(p^) (right), the 
lattice data and the corresponding fits (solid lines) to Eq. (f5) using the results reported in Tab. 2. 

The renormalized gluon propagator is then fitted to this last equation up to momenta p < 
4.638 GeV and the results can be found in Tab. 2 and seen in Fig. 1, where the lattice data 
and fitted curves appear displayed. We call the reader's attention to the good agreement for the 
fitted parameters from the two lattice volumes. This should be taken as a strong indication that 
finite-size lattice artifacts appear to play a not very important role in our present analysis. In 
addition, in the left plot of Fig. 2 we display the variation of the fitting parameters m?, and 
as a function of the fitting range, i.e. obtained from fitting the data up to a maximum value 
of the momentum, Pmax- We clearly observe a stability in the parameters' value from p ~ 2 GeV 
on. 
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(GeV2) 


(GeV2) 


(GeV^) 


(MeV) 




80 


0.7838(17) 


4.303(50) 


0.526(14) 


0.4929(39) 


338(3) 


1.31 


64 


0.7800(29) 


4.442(86) 


0.576(24) 


0.4964(64) 


334(4) 


1.72 



Table 2: Fitted parameters for Eq. (15) to the lattice data with the setup as in Tab. 1 and the frozen 
RGZ gluon mass, mo, defined by Eq. (7). 



2.3 A few words about the (A^) condensate 

Let us stress again here that the fit (15) is not exactly a tree level RGZ prediction, since Z ^ 1. 
Indeed, at tree level perturbation theory, we should have in principle Z^'^^ = 1, viz. the first 
order term of Z^'^^ = I + . . . + . . . g'^ + . . .. We could of course drop the Z by reestablishing 
the appropriate MOM subtraction scale via solving 

t^' + M^ _ 1 

with the fitted values of Tab. 2. Since we are fitting quite deep into the UV here, the extracted 
value for {A^) (or any other condensate) might emulate some of the already missing log 
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P„„„ [GeV] p [GeV] 



Figure 2: (Left) The value of fitting parameters obtained for p < p„„x, plotted in terms of p„,ax- (Right) 
the RGZ mass function given by Eq. (7) for the parameters of Tab. 2. The dotted lines account for the 
statistical errors to define the 1-sigma allowed region. The "tachyonic" contribution, resulting from the 
gluon condensate, is clearly dominating the mass function above roughly 0.5 GeV. 



information and as a consequence, become "unnaturally" large. Using the numbers of Tab. 2 
and the connection (3) we find 

(52^2),oz ^ 7 GeV2 (17) 

This estimate can be compared with the tree level OPE value of [51] that, obtained at a typical 
deep UV scale of 10 GeV, gives ~ 5(1) GeV^. The two aforementioned condensate estimates 
imply a serious violation of the roughly derived rule (12), suggestive of important missing in- 
formation. In particular, as we applied a tree-level RGZ formula to describe the lattice data 
and extract the condensate value, its natural renormalization scale should be the one obtained 
by solving Eq. (16), ~ 0.4 GeV, that lies deeply in the IR and makes questionable to use the 
renormalization group equation (RGE) to run the subtraction point from the OPE scale down 
to this momentum. Let us furthermore remark here that including RGE logarithms into Wilson 
coefficients during the OPE analysis is also of importance for a reliable determination of {A'^)^^^ 
(see, for instance, refs. [49,63]). 

For the rest of this paper, we shall not need to worry anymore about the subtleties of 
interpreting the fitting parameters in terms of condensates. We can simply take the high quality 
fit (15) for p < 4.638 GeV, see Tab. 2, and feed this into the GPDSE. Furthermore, as no 
important finite-size effect have been observed (the fitted parameters are compatible once their 
statistical errors are considered), we shall apply in the following the fit corresponding to data 
simulated at the largest lattice volume. 
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3 The ghost propagator DSE, its 1st iteration (one-loop pertur- 
bation theory), full numerical solution and the capital issue 
of the ghost-gluon vertex 

The goal motivating this section is the computation of the ghost dressing function by solving 
the ghost gap equation, 

where a{k'^) is the ghost self-energy that reads in d dimensions 

<k)=g J J^Hiik - q, k) _ D{q ) \^S^, ' (1^) 

the form factor Hi{k — q, k) is non- longitudinal one in the ghost-gluon vertex, 

Vf'iq - k, k; -q) = igor'^qu'Tu'uiq - K k; -q) 

= igor^mk-q),H^{k-q,k)-q,H2{k-q,k)] , (20) 

k — q and k are respectively the outgoing and incoming ghost momenta and go is the bare 
coupling constant. Now, the RGZ gluon propagator discussed in the previous section will be 
taken for D{q^) and applied to solve Eq. (18) either analytically, by taking F{q^) = 1 under 
the integral (as a first perturbative iteration), or by iterating numerically until a satisfactory 
convergence of results is reached. 

3.1 Analytical approach to the ghost dressing function 

We will take here F{q^) = 1 in the ghost self-energy of Eq. (19) and rewrite the RGZ gluon 
propagator obtained in the previous sections as 

^, 2x p'^ + M'^ ( a+ a_ \ 

+ (M^ -I- m'')p^ + A* \p^ + oj^ p-^ +uj±J 

where the poles are assumed to be cc, with positive real part. It has already been illustrated 
in [53] (SU(2) case) and [20, 52] (SU(3) case) that these forms allow for a consistent fit to large 
volume lattice data. Propagators with pairs of cc poles were already considered earlier in a DSE 
approach, see [64]. 

Now, on the ground of the OPE analysis performed in ref. [21], we model the non- longitudinal 
form factor, Hi{q, k) with the simplified non-perturbative formula: 

2 2 

H,iq,k) ^H,iq,0) =Wiq^) = 1 + , (22) 

based, as will be discussed below, on the approximation of a very soft dependence on the incoming 
ghost momentum, where 6 is a mass dimension scale discussed in ref. [21] and 

a' = ^(4^ 9'{A\,. . (23) 
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We expect that, even in the approximation of a soft dependence on the incoming momentum, 
the simple model vertex (22) can capture the non-perturbative essentials of the ghost-gluon 
vertex, while still allowing for analytical computations. In the numerical analysis performed in 
the next subsection, we will discuss a more general formula that has been derived in ref. [21] to 
describe more properly the — albeit not so precise — current lattice data [65]. Here we follow an 
approach, developed in [66] to study one-loop RGZ predictions for the ghost propagator, and 
as such we only summarize the underlying computation. The result with the tree-level vertex 
{W{q'^) = 1) has been previously obtained in the one-loop computation in ref. [67]. The one-loop 
contribution to the ghost self-energy coming from the model vertex (22) yields 

^/|f.W(*-^)Vl)^i>tf)(v-f^) , (24) 

where the non-perturbative extra piece for the vertex (that, being UV suppressed, will not induce 
extra UV divergences) is included under the integral. Writing 

0* + q^ + zo^ q'^ — ib^ 2b^ 
one eventually finds after some algebra for Eq. (19) the expression [66] 

(7(fc2) = z^g'^Nc [a+f{k',ul) + a-f{k\ojl) (26) 

+ a_/3+/(fc^wi^i62) + a_/?_I(fc^wi,-i62)j| , 

where ^ 

I{k\w\u^) = ^ - ..J,,,., (27) 



327r2e 1287r2A;4a;2 

2k^uj^{bk^ - 2z/2 + a;2) - (w^ _ 3^,2^4 ^ 3^2^4 ^ ^^i^2 _ 3^4^2^ 



-{k" + v^f ( In + In ^^^^ + 6kW In ^ 

+2V-fc4 - (1^2 _ ^2)2 _ 2A;2(jy2 + ^2) ((^^2 ^ ^2)2 ^ 2a;2(A;2 - + 

k^ -u^ + up- 



(ArcTan 
sj-k^ - (1/2 - a;2)2 _ 2fe2 {yi + Ji) 

k'^ + - u? 



-FArcTan 



- (zy2 _ ^2)2 _ 2p (^,2 + ^2)^ 



^It is assumed that Re [cj-'l > 0, next to i/^ e : 
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and we recall from the one-loop computation, with the tree-level vertex, in ref. [67] , the function: 



+ fcV(5fc2 + a;2 + A;2ln^-4A;2ln^^^)l ; (28) 

where everything has been computed in the MS scheme at the subtraction point ^. We now have 
all the ingredients to derive the ghost dressing function in the usual MOM scheme in d = 4. To 
this purpose, we need to properly renormalize the gluon propagator in the MOM scheme^. 
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D''''''{pX.^,. = J,. (29) 
and we simultaneously impose 

i^^°^(/)|,.=,. = 1 , (30) 

that can be achieved by the following (one-loop) subtraction again at the renormalization mo- 
mentum )U: 

where the ghost sclf-cncrgy, (j{k'^), in r.h.s. is given by Eq. (26) with zq determined by the MOM 
condition (29) and after replacing the coupling g by the one renormalized in the so-called Taylor 
scheme (see for instance ref. [51]). We will come back to the latter in the next subsection. 

3.2 Numerical approach to the ghost dressing function 

This section is devoted to deal with the computation of the ghost dressing function by the 
numerical resolution of the gap equation, Eq. (18), through a convergent iterative procedure. 
To this purpose, the bare ghost self-energy can be rewritten as 

.(,^) = ^oN. J - ,n (l - ^^) . (32) 

after the appropriate kinematical change {q ^ k — q) for the integration variable and particu- 
larization for d = 4. Then, after applying the MOM renormalization prescription to both ghost 
and gluon propagators, as in Eq. (31), the ghost gap equation reads 

MOM 



with 

■,2 



/ n o \ MUM 

[a{e) - ai^^')) (33) 



{a{k') - a(/x^))^°^ = ZlZ.^N, J q'dq K{q, k) F^°^{q') (34) 



^This condition implies a particular requirement for the overall factor, zo, in Eq. (21). 
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and the kernel given by 



where the subtraction procedure, as detailed in [24], is applied for two external momenta, k and 
p, chosen parallel and such that = fi^ , and Z^iZs) is the gluon (ghost) MOM renormalization 
constant. It should be noted that Hi in Eq. (35) is a bare but finite [68] quantity which needs 
no renormalization while, in front of the integral, the renormalization constants and the bare 
coupling especially appear in the right combination to cancel the UV singularities and give the 
MOM T-scheme coupling (e.g. see [51]), 

arif^') = ^Zi{f,')Z,{,,') . (36) 

This couphng aT(^t^) in Yang-Mills QCD is very well known from lattice data (e.g. see [51]). 
Thus, we will solve Eq. (33) by plugging the gluon propagator, modeled quite precisely with RGZ 
and lattice QCD as dccribed in sec. 2, into Eq. (35) to fully determine the kernel for the ghost 
DSE. Thus, we can solve the ghost DSE in isolation (i.e., without coupling it to the much more 
complicated gluon DSE) with the ghost-gluon form factor Hi as the only unknown quantity. 

In the previous analytical study, we plugged into the ghost self-energy a simple ansatz for the 
ghost-gluon form factor grounded on the OPE analysis of ref. [21]. Here, we will also consider 
the more general formula, 

Hi{q,k) Hi{q,0) ^ W{q') = c (l + ^J^) + (1 - c) ; (37) 

also grounded on the same OPE results, where estimates for b and defined by Eq. (23) can 
be borrowed from [21,51]. As explained in [21], c is a dimcnsionlcss parameter related to the 
perturbative contribution to the form factor Hi, which is usually assumed to be constant. The 
second term in Eq. (37)'s r.h.s. is included to keep the condition Hi{0,0) = 1, which appears 
to emerge, at least approximately, from 4D SU(2) lattice data [65]. Indeed, Eq. (37) is seen 
to capture the main features of the current (rough) lattice data for the ghost-gluon vertex. It 
should be also noticed that the simpler form factor used in the previous subsection, Eq. (22), 
will be recovered from Eq. (37) as c = 1. 

Thus, Eq. (22) and Eq. (37) can be applied in Eq. (32) to solve the gap equation (33) with 
the lattice gluon input, and the parameters c and w taken as free parameters that must be sized 
to properly account for ghost dressing lattice data. This will be done in next subsection. 



3.3 Results and compcirison 

Let us now discuss the solutions of the GPDSE for the two classes of ghost-gluon vertices (22) and 
(37) considered here. Recall that both vertices come from the OPE analysis of the corresponding 
pure Yang-Mills vertex [21,51]. As will be shown, the first vertex, being an approximation to 
the second vertex by neglecting the perturbative contribution at the subtraction point (c = 1), 
implies that the OPE condensate (included in the parameter o^) to mimic this contribution, 
resulting in an enhanced value for a^. Our results are summarized in Fig. 3 and Tab. 3. 
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In the following, when computing the ghost self-energy, we take axip = 3.6) = 0.25 from 
ref. [51]. Furthermore, the different outcomes of the GPDSE are confronted to the ghost dressing 
lattice data published in ref. [36], simulated at 80^ and 64^ lattices at P = 5.7. 

As a self-consistency check, we have found it worthwhile to display also in each case the vertex 
used, to facilitate a comparison with, for instance, the SU(3) vertex lattice data of ref. [69] or 
future more precise results. 

3.3.1 Case 1: a trivial vertex 

Our first results are those of the fully iterated GPDSE for the usually made trivial approximation 
Hi = 1 (named HIO) and in addition, we also consider a "boosted by hand" constant behaviour 
for the form factor (named Hlctc); see in Fig. 3 and Tab. 3. We observe that, as was pointed 
in refs. [15, 24], neither approximations are able to properly account for the behaviour of lattice 
data. Indeed, the tree level vertex HIO shows deviations for momenta below 2 GeV, while 
the Hlcte ghost dressing function, although being closer to the lattice data, also shows clear 
deviations from the corresponding lattice function. 

For comparison, we have also shown in Fig. 3 the once iterated approximation which, as 
we have explicitly checked, corresponds one-loop perturbation theory using the RGZ formalism, 
with the same vertices HIO and Hlcte (labeled as HlOliter and Hlcteliter). Again, deviations 
form the lattice ghost data in the infrared region are striking. 

3.3.2 Case 2: RGZ and NP solutions with vertex model (22) 

As a next approximation, we take into account the simple vertex model (22) (or equivalently 
the vertex model (37) with c = 1) that can describe reasonably well the vertex lattice data. We 
will reconsider both the single iteration (named RGZ) and non-perturbative (NP) results of the 
GPDSE to determine, in each case, the optimal a and b to the ghost propagator data by treating 
them as fitting parameters. As a cross-check, for the single interation results, we applied both the 
numerical integration and the closed form expression (26) and found no significant discrepancy 
(less than 0.5 % and only for very low momenta, < 0.03 GeV). The results can be found in 
Tab. 3. Again for the sake of comparison, we also plotted in Fig. 3 the results for the single 
iteration procedure with optimal a and b corresponding with NP results (labeled as NPlitcr). 

The vertex (22) with the derived value of a and b is shown in Fig. 4. It is instructive, using 
the optimal a of both cases, to see what value of {A^)^^^^ comes out, something that can be 
achieved using the OPE motivated formula (23). We get (^^)ope - 22.6 GeV^ for NP and 
{A'^)qp^ ~ 23.9 GeV^ for RGZ which are both unnaturally large when compared to the result 
of the OPE analysis of Taylor coupling in ref. [51]. This is a manifest consequence of the neglec- 
tion of the perturbative contribution to the ghost-gluon vertex at the subtraction point (choice 
c = 1) which is effectively borrowed by the OPE condensate. It should also be noticed that, for 
obtaining a reasonably proper description of lattice data with the single iteration results, the 
vertex needs to get quite strongly enhanced at low momenta, mainly due to the small size of the 
parameter b. 

To get an idea what the role is of fully iterating the GPDSE compared to the single iteration 
of the previous case, we have also shown the once iterated version using the same vertex (and 
the same parameters) as for the NP case (labeled with NP liter). Manifestly large discrepancies 
start to appear for momenta below around 2 GeV and imply that the higher order corrections 
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to the one-loop RGZ ghost propagator are relevant when precise results are the objective. The 
same could be concluded by comparing HIO and HlOliter cases in the left plot of Fig. 3 and 
Hlcte and Hlcteliter in the right plot (solid and dotted blue lines). 

3.3.3 Case 3: fully iterated GPDSE with vertex model (37) 

With the full solution of the GPDSE using the somewhat more general vertex model Eq. (37), 
our primary goal is to present both a better description of lattice data and relate as much as 
possible the vertex parameters to independent analysis in order to give them some "physical" 
meaning. Then, as the full iteration makes it rather difficult to treat all parameters as to be 
fitted, we opt to borrow a from the relation Eq. (23) and we take g^{A^) ope — 7 GeV^, at tree- 
level and our renormalization point 3.6 GeV, from ref. [51]. We then look for the optimal values 
of c and w, by the minimization of x^/d.o.f. when fitting the ghost DSE prediction to the ghost 
dressing lattice data, for several values of b around the SU(2) one reported in ref. [21] and then 
we simply took the best b. We notice that changing b corresponds to either shifting the vertex 
peak towards the large-momentum region and reducing its height or shifting it towards the low- 
momentum and increasing its height, thence provided that o^, i.e. {A'^), is fixed, the value 
of b should not change too much. Fig. 3 illustrates clearly that this GPDSE solution (labeled 
HIOPE) captures the lattice data over the full momentum range very well, with the parameters 
shown in Tab. 3. Also for the vertex itself, see Fig. 4, we observe a good resemblance^ with the 
SU(3) lattice vertex of ref. [69] (see Fig. 4 of this paper) or with the SU(2) results of [70]. 



label 


c 


a 


b 


w 


xVd.o.f. 






(GeV) 


(GcV) 


(GeV) 




HIOPE 


1.26 


0.80 


1.3 


0.65 


0.83 


Hlcte 


1.35 









13.0 


NP 


1 


1.46 


1.4 




1.74 


RGZ 


1 


1.50 


0.85 




1.98 


HIO 


1 












Table 3: Parameters for the form factor Hi, given by Eq. (37), used to solve the ghost DSE. The set of 
parameters are refereed in Fig. 3 as labeled in the first column. The first three lines correspond to fits 
over the whole momenta window, while the fourth stands only for momenta above 0.3 GeV. 



4 Conclusion 

As the main goal of this paper, we have shown the quantitative compatibility between large- 
volume lattice ghost propagator and its DSE description. To this purpose, we discussed the 
relevant role played by a non-trivial ghost-gluon form factor (in particular, the non-longitudinal 

^It should be noticed that the authors of ref. [69] plotted Z^^, which corresponds to the bare but finite form 
factor Hi , rcscalcd to bo 1 at 3 GoV (and dropped any possible large-momentum purely pcrturbativc contribution 
away, as the one coded by c in Eq. (37)). Although, for a proper comparison with our vertex results, we would 
have required the non-rescaled bare form factor, but one can anyhow realize that the kinematical structure of 
their form factor matches pretty well with our result. 
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p [GeV] p [GeV] 



Figure 3: (Left) ghost dressing data obtained by solving tlie ghost DSE, Eq. (33), with a ghost self-energy 
given by Eq. (34) and the gluon propagator predicted from Eq. (5) and fitted to lattice data as input for 
the kernel in Eq. (35). The results for the different set of parameters for the form factor Hi collected 
(and labeled) in Tab. 3 are confronted to the ghost dressing lattice data taken from ref. [36]. (Right) the 
same but applying the analytical perturbative approach described in subsec. 3.1. The meaning of the 
labels for the different curves is explained in the main text. 

form factor usually called Hi) that has been recently modeled founded in the OPE vertex de- 
scription and in good agreement with current (although not yet very precise) vertex lattice 
results. Indeed, we have shown that, only after the appropriate inclusion of such a non-trivial 
form factor, the integration of the gluon propagator in the ghost DSE will result in a ghost 
propagator numerically consistent with the large-volume lattice data. A key point in our anal- 
ysis is the fact that, once the gluon propagator is precisely borrowed from the lattice and the 
renormalized coupling unambiguously related to the one in Taylor scheme, the only unknown 
DSE ingredient is the above form factor. 

As an easy-to-handle model for the gluon propagator is of a great help to solve the ghost DSE, 
after being properly fitted to lattice data, we invoked the tree-level gluon propagator within the 
RGZ approach. This RGZ gluon propagator fulfills a two-sided goal: it provides, after applying 
the usual MOM renormalization prescription, with an appropriate model describing lattice data 
to be plugged into the DSE and solve it numerically; and, being so simple as a model, it allows 
for an approximate analytical result based on the one-loop RGZ approach. Then, the results are 
compared and we concluded that the perturbative RGZ approach (= 1st iteration of GPDSE) 
allows for a good description of the non-perturbative lattice results only after one considerably 
enhances the ghost-gluon vertex at low momenta, a fact not really supported by the lattice data 
of [69]. 
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Figure 4: The different vertex models used to obtain the numerical or the single iteration results for the 
ghost dressing function. They correspond to Eq. (37) where the different set of parameters, labeled as in 
Fig. 3, should be read off from Tab. 3. 

and Adelaide lattice groups for providing us with their data. 
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